Casimir repulsion between Topological Insulators in the diluted regime 
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The Pairwise Summation Approximation (PSA) of Casimir energy is applied to a system of two dielectrics 
witii magnetoelectric coupling. In particular, the case of Topological Insulators (TI) is studied in detail. De- 
pending on the the optical response of the TI, we obtain a stable equilibrium distance, atraction for all distances, 
or repulsion for all distances at zero temperature. This equilibrium distance disappears in the high temperature 
limit. These results are independent on the geometry of the TI, but are only valid in the diluted approximation. 
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I. INTRODUCTION 

Since 1948, when Casimir introduced the energy that got 
his name [I], calculation formulas have been looked for 
Many analytical and numerical methods have been proposed, 
such as the zeta function technique, the heat kernel method, 
semiclassical methods or Green function (local) methods just 
to mention a few of them 1 2 1 . However, exact results have 
been obtained only for some simple geometries. In this con- 
text, several approximations have been proposed such as the 
Proximity Force Approximation or PSA, to cite two of them. 

More recently, a novel Multiscattering formalism of 
Casimir effect for the electromagnetic (EM) field have been 
proposed 13] |4j. It had been successfully used in calcula- 
tions of Casimir effect in new geometries, but has also been 
applied to obtain a generalization of PSA to dielectrics with 
magnetic response in |5 1. In addition to that, this Multiscatter- 
ing formalism was applied to obtain a generalization of Earn- 
shaw's theorem to Casimir effect, claiming on the impossi- 
bility of equilibrium states in systems governed by Casimir 
physics |6|. 

On the other hand, it has been suggested that TI could lead 
to the appearance of equilibrium systems because Casimir ef- 
fect in [71 . 

Our goal is the generalization of PSA to obtain the Casimir 
energy between dielectrics with magnetoelectric response. To 
do so, we use the generalized constitutive relations of materi- 
als with magnetoelectric coupling to obtain a potential, which 
is the tree level of the Born series of the required T matrix. 

We will apply this formalism to the Casimir energy between 
TI, as a result, we obtain a prelora of different behavior of the 
system depending on the electromagnetic response of the TI. 
Atraction for all distances, repulsion for all distances, or also 
the existence of a stable equilibrium distance are obtainable in 
this system at zero and finite temperature. 

Then we find a system where the recent extension of Earn- 
shaw's theorem of inexistence of stability on systems gov- 
erned by Casimir physics |6| is not applicable (magnetoelec- 
tric coupling were not considered in |6|), therefore stable po- 
sitions are not forbidden in the system studied here. 

In addition to that, we also propose a system candidate to 
perform quantum levitation in vacuum. Other examples of re- 
pulsion within Casimir systems have been proposed for very 
specific geometries in vacuum f8l|, for systems immersed in 
fluids with dielectric properties different of vacuum |9 1, or for 
the Casimir interaction between dielectrics and metamateri- 
alslHOj. 



In this article we extend the results for TI for arbitrary ge- 
ometry in the PSA, valid only in the diluted limit. It is, when 
the dielectric constant e^ — > eo and the diamagnetic constant 
Ha —>■ Mo for ^ny ct body of the system. 

The article is structured as follows; Using the generalized 
constitutive relations of materials with magnetoelectric cou- 
pling and the formalism presented in f5|, in Sect. |ll]we will 
obtain the PSA of Casimir energy for the zero and high tem- 
perature limits starting from the exact Casimir energy formula 
given in f3l. In Sect. [Ill] we will apply the formalism devel- 
oped in Sect. Illlto obtain the PSA Casimir energy for TI in the 
zero and high temperature limits. We will obtain conditions to 
the appearance of atraction, repulsion and the appearance of 
equilibrium distances in term of the response of the TI to the 
EM field. Finally, in Sect. IV we will perform a numeri- 



cal study of the PSA Casimir energy between two TI parallel 
plates to compare the results of PSA with exact results already 
done in f7|, to study the validity of the approximations made. 
We also will obtain numerical results of the PSA Casimir en- 
ergy for the sphere-plate system, because its experimental rel- 
evance ITTl . 



II. GENERALIZATION OF PSA TO DIELECTRICS WITH 
MAGNETOELECTRIC COUPLING 

Our goal is the calculation of the complete electromagnetic 
Casimir energy between two bodies in the soft dielectric limit 
with magnetoelectric response. For this purpose, we use the 
multiscattering approach to EM Casimir effect |3 , 4| to obtain 
the Casimir energy between two compact bodies at a given 
temperature T as 



Et ^ ksT^ ln|l -N(k„ 



(1) 



n=0 



are the Matsubara 



where 1 is the identity matrix, k„ — 

frequencies and At — 2^1^ t ^^ ^^ thermal wavelength. The 
prime indicates that the zero Matsubara frequency contribu- 
tion has height of 1/2. All the information regarding the sys- 
tem is described by the N matrix. For a system of two objects, 
this matrix is N = T1U12T2U21. T^ is the T scattering ma- 
trix of the rth object, which accounts for all the geometrical 
information and electromagnetic properties of the object. Uij 
is the translation matrix of electromagnetic waves from ob- 
ject i to object j, which accounts for all information regarding 
the relative positions between the objects of the system. From 



Eq. ([T]), the quantum (T -^ 0) and classical limits (h -^ 0, 
equivalent to the high temperature limit) are easily obtained 
as 



En^ 



he 

2^ 



dfcln|l-N(K)| 



and 



ErI — 



kuT 



ln|l-N(0)| 



(2) 



(3) 



respectively. To obtain the PSA of Casimir energy, we ap- 
ply the formalism developed in IJS). Then, using the relation 

\iy\A\ = Tr(ln(yl)) and that ln(l - x) = - E^i ^' we 
write Eq. ([TJ as the asymptotic series 



CO ' CO _. 

i?T = -fcsrV V-Tr(NPK)) 



(4) 



n— p— 1 



In addition to that, we use the position representation of op- 
erators instead the multipole representation used in |[3), so we 
can identify U^^ = Gqq^, where Gq^p is the matricial free 
dyadic Green function given in |5|. 

The T operator is related with the potential V by the 
Lippmann-Schwinger equation. When V is small, we can ap- 
ply a Born expansion to the Lippmann-Schwinger equation 
to obtain an approximation for the T operator in the diluted 
limit as T^ sa Vi. Here we remind that the use of a Born ap- 
proximation as a result of the Lippmann-Schwinger equation 
is more valid for lower potentials. In our case this means that 
PSA will be valid in the diluted limit, and we will need high 
order corrections when this diluted limit does not be longer 
applicable. 

To define the potential 1/,, we use the generalized constitu- 
tive relations of materials with magnetoelectric coupling 



D = eE + aH, 

B = /3E + iM, 



(5) 



where a and /3 are the magnetoelectric couplings. The poten- 
tial of each body can be defined as the difference of energy of 
the EM field because the existence of this body. Having into 
account that the energy of the EM field is defined as 



E=]: I dx^ (E • D + H • B) , 

J \l 



(6) 



where D = eqE and H = jiqR in the vacuum, the EM energy 
in presence of N generalized dielectrics is 



1 /• ^ 

E=- dx'' (eoE^ + ^oH") + Yl ^^- 
•^^^ 1=1 



(7) 



We use the generalized constitutive relations given in Eq. Q 
to obtain the excess of energy because the existence of each 
dielectric as 



AE, 



da;^(E,H) 




then the potential is defined as 



V, 



fJ^i 



Xi (r) , 



(9) 



where li = e^— eg, jli — Hi—Ho and Xi (r) is the characteristic 
function of the i body volume (1 inside the body and in 
the rest of the space). The generalization to space dependent 
dielectric constants is straightforward. 

Now we study the lowest expansion order in dielectric 
quantities of the Casimir energy. In the lowest expansion or- 
der (pMax = 1 and T,; — Vi in Eq. Q), we get the asymptotic 
approximation of the Casimir energy between two bodies as 



El 



OO 



~fcBr^Tr(l/iGo42l^2Go,2l): 



(10) 



where all operators depends on k„. These operators are de- 
fined over three different linear spaces fS): (1) An EH space, 
whose components are the electric and the magnetic field; (2) 
over the space coordinates, because we are working with a 
vector and a pseudovector; and (3) over positions. We must 
solve the trace over these three spaces: EH space, vector coor- 
dinate space, and position space and sum over the Matsubara 
frequencies to obtain the PSA energy. 

As a first application, we use the potential Vi given in Eq. 
(|9]l to obtain the PSA Casimir energy of frequency indepen- 
dent dielectric with magnetoelectric couplings at zero temper- 
ature as 



E = 



-he 



{A^y 



■ la 



dridr2 
1J2 |ri-r2r 



(11) 



where 



7o = 23eie2 - 7iiil2 ~ Ijiih + 23/ii/i2 
+ 7aia2 + 23ai^2 + 23;3ia2 + iPih, 



(12) 



This result generalizes the Feinberg and Sucher potential lfT2l 
to objects with magnetoelectric couplings, but only in the di- 
luted limit. In the high temperature limit, the PSA of Casimir 
energy is 



E^ 



-kuT 



(47r)- 



-Id 



dridY2 

2 |ri-r2r 



where 



7ci = 3eie2 + 3/2i/i2 + 3ai^2 + 3/3ia2. 



(13) 



(14) 



It is a straightforward calculation obtain the PSA at any fi- 
nite temperature. The result is similar to the finite temperature 
results of PSA for diluted dielectrics shown in [5|. 

As we can see, if magnetoelectric couplings have different 
signs and are strong enough to compensate the usual electric 
and magnetic coupling of dielectrics, a repulsive Casimir en- 
ergy between objects can be achieved in the diluted limit. 



III. PSA FOR TI 

Recently, in Q, A. Grushin and A. Cortijo demonstrated 
the existence of an equilibrium distance between topological 
insulators with opposite topological polarizabilities sign be- 
tween parallel plates. This result was extended for all tem- 
peratures in 1 13 1, leading to a reduction of the equilibrium 
distance with the temperature until the disappearance of the 
equilibrium distance in the high temperature regime. Regions 
of Casimir repulsion for all distances were also reported for 
some values of w = -^^ and \9\. In this part of the article, 
we apply the formalism developed in the previous section to 
TI. The constitutive relations of TI are D = eE + a {d/n) B 
and H = B//i — a (O/tt) E. In order to apply the formalism 
developed in the last section, we need the fields D and B as a 
function of E and H, then 



D= e 



B 




/la [ - ) E + fj,H. 



(15) 



where Ui = a—, w = ■^^, x = ^^^Iri — rol is the dimen- 

i TT ' LOR ' C I -"- ^1 

sionless distance and fi{x) and f2{x) are functions which 
depends on a linear combination of different Meijer G func- 
tions too long to be written here. All terms of Eq. ([T9| are 
positive ifi{x) and f2{x) are smooth positive functions in 
the real positive axis, see Fig. [TJ for any value of di and 
62, except the product aia2, which becomes negative when 
sign(6'i) = -sign(e2)- 

At short distances we can approximate /i(x) sa Stt and 
/2 (x) ~ 6tt, then 70 can be approximated by 



7o = 60aia2 + 1Zoi\a^. 



(20) 



If sign(6'i) — — sign(02), the energy Eq^^ becomes positive 
(i.e. repulsive) for all 



'h02< 



607r2 
23^ 



(21) 



For diluted materials, this condition is beyond the positive en- 
ergy condition 1^1021 < ^^£1(0)62(0) < i^ na, then 
we can consider it is always fulfilled for diluted dielectrics. 



where a is the fine structure constant (a = fj) Following the 
discussion in ifTJl , we assume fi ^ jiq and 9 is quantized in 
odd integer values of it such that 



= (2n+l)7r, 



(16) 



where n e Z, determined by the nature of the magnetic coat- 
ing, but independent of the absolute value of the magnetiza- 
tion of the coating. Positive or negative values of 9 are related 
to different signs of the magnetization on the surface I.14J . 
which we consider is perpendicular to the surface of the body. 
Being a topological contribution, 9 is defined in the bulk as 
a constant whenever the bulk Brillouin zone is defined ifTSl . 
In addition to that, we assume that the frequency dependent 
dielectric function e{uj) is described by an oscillator model of 
the form: 



e{iK) = eo + 






'r.i + 7fl,iCK + C^K^ 



(17) 



Because there are few experimental results of e{iK) for TI, 
we will assume that there is just one resonance. The case of 
multiple resonances or another model of e{iK) can be easily 
generalized [7|. 

Because the frequency dependence of eiin) for TI, it is not 
possible to obtain the integrand of the PSA energy at zero tem- 
perature in a closed form in terms of simple functions. Then 
the PSA Casimir energy at zero temperature for topological 
insulators is given by the integral 

he f f dridr2 ,.p, 

f7o, (18) 



j^PSA 
Eq 



(47r)3 



r2| 



where 70 now depends on the distance between points of the 
TI and on the dielectric properties of both TI as 



7o — w\w2xfi{x) + 60Q;ia2 + ^ia^a^ 



+ {a\w1 + a^wl) xf2{x), 



(19) 




Figure 1. (Color online) fi{x) (red curve) and f2{x) (black curve) 
as a function of the dimensionless distance x. The blue curve is their 
common asymptotic limit to — at large x. The limits at small x of 
fi{x — >■ 0) — ^ Stt (yellow curve) and f2{x — >■ 0) — ^ Gtt (green 
curve) are also plotted. 



At large distances we can approximate fi{x) w — and 
f2{x) w ^, then 7o can be approximated by 

O l^ Q Q O Q Q O \J\J 

7o = W1W2 + a^a2 + aiW2 + 0^2""^! + ^"i«2- (22) 



If sign(0i) — — sign(02) and wi = W2 
comes positive (i.e. repulsive) for all 



w, the energy be- 



2 23 

W < ITT. 

46 



, , 240 46 _„ „\ 

+ "2 ^"i"2 - ^"l«2 I 

(23) 

As a consequence, there are several different regimes for 

the system in the quantum limit. When sign(6'i) = sign(02)> 

the Casimir energy is enlarged because the contribution of 



topological polarizability 9, but when sign(6'i) = — sign(6'2), 
different regimes appear 

When sign(6'i) = — sign(6'2) and for low enough absolute 
values of Oi and 62, the condition given in Eq. ( |2T] i is fulfilled 
while the condition given in Eq. ( |23| l is not, then we have a 
repulsive Casimir energy at short distances and an attractive 
Casimir energy at large distances, then there must exist a sta- 
ble equilibrium distance in this case. 

But if the condition given in Eq. (|23| is also fulfilled, 
the magnitude of the positive arguments of Eq. ([T9| is not 
large enough to compensate the repulsion because topological 
charges 6, then we obtain repulsion for all distances. 

In the high temperature limit, another result is obtained. In 
this case the PSA gives the Casimir energy as 



E. 



PSA 



ikeT 



dridr2 

2 |ri-r2r 



(24) 



where, contrary to the zero temperature case, -fci does not de- 
pend on the distance between points, and it is given by 



7ci 



2„,,2 



W^ W. 



— 2 2 —2 2 

aiW2 + Ct2'^\ 



a](X2 + 2Q;ia2. (25) 



Depending on the values of w and j^j, we obtain or repulsion 
for all distances or atraction for all distances. The condition to 
be fulfilled in order to obtain repulsion in the classical limit, 
when w\ — W2 — w v& 



w^ < 



1 



+ ttj — 8aia2 — 2a\ai 



(26) 
In Fig. [2] we represent the behavior of the PSA energy as a 
function of 6 and w when wi = W2 and di = 9 = —62 in the 
quantum and in the classical limit, having into account there 
exist a forbidden region because the positive energy condition 
fTFl. In both cases we find different regions of parameters of 
repulsion for all distances, but in the quantum limit there is 
a region of existence of an equilibrium distance, while in the 
classical limit a region of atraction for all distances appears. 

It is expected that, at finite temperatures, the equilibrium 
distances found in the quantum limit would be reduced when 
temperature increases |13|, and Fig. l2] suggest that behavior, 
but we have let this study as a future work. 



IV. CASIMIR ENERGY BETWEEN SPHERES AND 
PLATES 

In this section we will use the PSA formulas presented in 
the previous one to obtain the Casimir energy to specific sys- 
tems. In particular, we study the sphere-plate system because 
it experimental relevance and the two infinite plates system to 
compare the PSA results with the exact ones. Both systems 
were studied at zero temperature. 

We were not able to obtain analytical results, so a numerical 
integration procedure has been implemented. 

The two infinite parallel TI plates has been already studied 
in Q. Our interest here is to compare the results obtained by 
the PSA with the exact ones. We impose 9i = — 6*2 = tt in 




200 



Figure 2. (Color online) Atraction versus Repulsion in the classical 
limit (T — >■ 00) and in the quantum limit (T — >■ 0) as a function 
of w and \6\. For the Classical limit, the green region shows repul- 
sion for all distances and the rest of phase space shows atraction for 
all distances. For the quantum limit, the red region shows repulsion 
for all distances and the white region shows parameters for which 
an equilibrium distance appears. Several curves of constant equilib- 
rium dimensionless distance x^q are also plotted. The black region 
is a forbidden regi on fo r the parameters due to the positive energy 
condition \a\ < ^/e(0) (16). 



order to obtain an equilibrium distance, and different w from 
w = 0.2 to w — 0.6. The obtained results of PSA and the 
exact ones are compared in Fig. [3] where we find that PSA 
tends to overestimate the magnitude of the Casimir energy the 
more the larger w. This fact reflects the nature of the approxi- 
mation made in PSA, where we must assume the diluted limit 
is valid. On the other side, an excellent approximation of the 
equilibrium distance between plates is obtained, then the PSA 
gives a good qualitative result. 

We have also used the PSA to study the sphere-plate sys- 
tem because it experimental relevance in Casimir effect ex- 
periments. In this case, we also impose 9i = —6*2 — tt and 
w — 0.45 for both plate and sphere in order to obtain an equi- 
librium distance at T = (see Fig. |2]). We vary the dimen- 
sionless radius of the sphere Rg — '^^Rs from ^^ — > until 
Rs = 1 to observe it effect in the equilibrium distance. In 
Fig. |4] the PSA energy per unit of volume of the sphere is 
plotted as a function of the distance between nearest points of 
the plate and the sphere. As a result, an equilibrium distance 
has been obtained for all studied sphere radii. This equilib- 
rium distance reduces when the radius increases until reaching 
a constant value. The appearance of an equilibrium distance is 
easy to understand because the nature of the integrand in PSA 
calculations at zero temperature. As discussed in Eq. ([T9|, 
from a given distance between sphere and plate, the points of 
the sphere nearer to the plate tends to increase the Casimir en- 
ergy (giving a repulsive contribution), while the rest of points 




Figure 3. (Color online) Comparison between exact Casimir energies 
in units of Eo = Ahc/{2Tv)^ {ojR/cf(M\ curves) and PSA (points) 
at zero temperature for TI infinity plates with 6i — —62 ~ n and 
different e(0) as a function of the dimensionless distance d 1 7 1. PSA 
tends to overestimate the absolute value of the energy, but captures 
the equilibrium distance without appreciable error. The overestima- 
tion of the energy increases with e(0), as expected because the nature 
of the approximation. 



of the sphere tends to reduce this energy (giving an atractive 
contribution). As at short distances Eq. ([T9| diverges as Jr 

and L ^ diverges, it is evident that there always be a distance 
where repulsion compensates atraction, reaching the system 
an equiUbrium distance. 
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Figure 4. (Color online) Dimensionless PSA Casimir energies in 
the sphere-plate TI system per unit of volume of the sphere e = 
-^ ■§ with 9i — —62 ~ TT and w = 0.45 as a function of the 

dimensionless distance between the plate and the contact point of the 
sphere h — d — Rs, where d is the dimensionless distance from the 
surface of the plate to the center of the sphere. Positive energies are 
the whole curves, while negative energies are the dashed curves. The 
equilibrium distance for each sphere radius is represented by a dot. 
These dots are joined by a soft line to have an idea of the dependence 
of the equilibrium distance with the radius of the sphere. 



V. CONCLUSIONS 

In this article we have extended the PSA to dielectrics with 
magnetoelectric coupling. The extension to anisotropic ma- 
terials of this formalism is straightforward, although we have 
not done it here, is shown how to do it for the case where it 
would be needed. 

There exist several examples of magnetoelectric materials, 
as metamaterials flOl, already proposed as candidates of re- 
pulsive Casimir effect, and TI 171, which have been proposed 
recently to obtain equilibrium systems in Casimir physics. 

We remind here that the extension of Eamshaw's theorem 
to Casimir effect |6| does not forbid the appareance of equi- 
librium distances between TI, because magnetoelectric cou- 
plings were not considered. 

As a result, we have demonstrated that two TI or arbitrary 
shape with opposed signs sign(6'i) = — sign((?2) can present 
or repulsion for all distances or the appearance of a stable 
equilibrium distance at zero temperature as a function of their 
dielectric properties, given by parameters wy \d\ in the model 
studied here. 

At high temperatures the behavior of the system differs. 
Equilibrium distance collapses to zero, and the repulsion pa- 
rameters range is reduced. Then we obtain, as a function of 
the dielectric properties of the TI, or atraction or repulsion for 
all distances. 

We also have used the results derived here to study the 
Casimir energy for a system of two infinity TI plates, to com- 
pare the PSA with exact results given in |7|. As a result, we 
found that PSA gives better quantitative results for lower w, 
as expected because PSA works better in the diluted limit. 
On the other side, we obtain an excellent approximation of 
the equilibrium distance between plates, then we obtain good 
qualitative results for all 9 and w studied. 

We have also used the PSA to study the Casimir energy be- 
tween an infinite TI plate and a TI sphere because the exper- 
imental relevance of this geometry in experiments of Casimir 
effect [11 J. 

The results presented in this article are valid for objects 
of arbitrary shape, but in the diluted limit. Even so, this is 
the relevant limit because topological polarizabilities 9 scale 
with the fine structure constant a, and Casimir repulsion and 
Casimir equilibrium are more relevant when e is small too, in 
the diluted limit itself. 
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